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Saturation effects affecting absorption and fluorescence spectra of an atomic vapor confined in 
an Extremely Thin Cell (cell thickness L < l/j,m) are investigated experimentally and theoretically. 
The study is performed on the D2 line (A = 852nm) of Cs and concentrates on the two situations 
L = A/2 and L = A, the most contrasted ones with respect to the length dependence of the 
coherent Dicke narrowing. For L — A/2, the Dicke-narrowed absorption profile simply broadens 
and saturates in amplitude when increasing the light intensity, while for L = A, sub-Doppler dips 
of reduced absorption at line-center appear on the broad absorption profile. For a fluorescence 
detection at L — A, saturation induces narrow dips, but only for hyperfine components undergoing 
a population loss through optical pumping. These experimental results are interpreted with the 
help of the various existing models, and are compared with numerical calculations based upon a 
two-level modelling that considers both a closed and an open system. 
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I. INTRODUCTION 

High resolution spectroscopy of atoms confined in a 
thin cell is promising for the investigation of complex 
spectra of atoms and molecules. When irradiating a 
10 — 1000/im thin cell under normal incidence, the 
transmission spectrum of a single light beam had re- 
vealed tiny sub-Doppler features, nevertheless allowing 
the resolution of the hyperfine structure (h.f.s.) of al- 
kali atoms resonances (JUL 13- Indeed, the atoms travel 
wall-to-wall and resonantly interact with light according 
to a transient regime, yielding an enhanced response for 
atoms with longer interaction time, i.e. for those atoms 
insensitive to a Doppler shift. Recently, the extension 
to a technology of Extremely Thin Cell (ETC, thickness 
L < 1/im) [J] has shown that these sub-Doppler reso- 
nances become easily observed and can be a dominant 
effect, even in the absence of a frequency modulation 
(FM) technique. With these ETCs, an even narrower 
fluorescence spectrum can also be detected, showing that 
transmission and fluorescence behave differently. For the 
transmission spectrum, it was shown that one observes a 
maximal coherent sub-Doppler Dicke narrowing [f| for a 
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cell length L = A/2, with periodic revivals at (2n + l)A/2 
0, 0] , while for L = A (or nX ) , the local coherent atomic 
response is phase- mismatched all over the ETC, mak- 
ing the sub-Doppler features vanishing and only allow- 
ing a residual Doppler-broadened response. This length- 
dependent narrowing has no equivalent in the detection of 
an incoherent process such as the fluorescence 0, whose 
narrow features simply originate from the enhancement 
of the contribution of those atoms whose velocity allows 
for a long interaction time. This leads to a behavior that 
is monotonic with the cell thickness. These features of 
sub-Doppler transmission and fluorescence in ETCs were 
demonstrated in the linear regime of interaction (with re- 
spect to intensity), allowing for example the experimental 
determination of transition probabilities of the respective 
h.f.s. components @|. 

The detailed mechanisms for these sub-Doppler fea- 
tures are known to be complex because, depending upon 
the irradiating intensity, a coherent linear regime [2| and 
an (incoherent) non-linear regime of optical pumping had 
already been distinguished for the absorption in relatively 
long cells (L > 10/xm)[]J,|2j,l3j]. As long as the atomic sys- 
tem can be described in a two-level frame, various treat- 
ments for transmission experiments were developed in an 
asymptotic regime relatively to the irradiating intensity 
[E 0- 0, [§] , up to a full analytical formal expansion ([icl|). 
For ETCs, these high-intensity effects become even more 
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complex, with respect to the interplay between the inter- 
ferometric dependence associated to the Dicke narrow- 
ing, and the velocity-dependent efficiency of saturation 
mechanisms. 

In this communication, we present an experimental 
and theoretical study of the intensity effects in an ETC. 
For simplicity, we restrict our experiments to the two 
cases L = A/2 and L = A, the most important ones with 
respect to the periodicity of the Dicke coherent narrow- 
ing. An additional simplicity for this restriction is that 
it allows to neglect some of the Fabry-Perot effects in- 
trinsic to ETCs, namely the mixing of transmission with 
reflection signals because the non resonant reflection 
vanishes for these cell lengths. The study is performed 
on the D2 line of Cs vapor, with a spectral resolution 
about an order of magnitude better than in , and with 
an irradiation intensity orders of magnitude larger than 
in previous experiments on Cs vapor in an ETC [|| 0] • 
In addition to the known occurrence of narrow dips over 
a broader background in the absorption spectrum 0, [3 j 
we discuss here the appearance of narrow dips over the 
(sub-Doppler) fluorescence spectrum. Such dips of re- 
duced fluorescence were only briefly described in prelim- 
inary reports [l2T |. Here we present experimental and 
theoretical studies aimed at the clarification of the ori- 
gin of the observed narrow dip in the fluorescence pro- 
file. An interesting peculiarity of the observed narrow 
reduced-fluorescence dip is that it appears for all h.f.s. 
components but the one that does not suffer population 
loss due to hyperfine and Zeeman optical pumping. 

In spite of the complexity of the Cs atomic system with 
respect to saturation effects (as due to the many hyper- 
fine and Zeeman substates), we show, numerically as well 
as on the basis of general theoretical arguments, that the 
major features of our observations can be interpreted in 
the frame of a two-level model provided that closed and 
open atomic systems are distinguished. 



is controlled by the temperature T of the side arm (un- 
less stated otherwise, T = 119°C). The irradiating beam 
has a diameter of 0.4mm. Its intensity is controlled with 
neutral density filters F±. The transmitted light power 
is measured by the photodiode PD1. To ensure a con- 
stant sensitivity of the detector, the off-resonance inten- 
sity falling onto PDl is kept constant by filters F2. To 
record fluorescence spectra, the photodiode PD2 collects 
the induced fluorescence emitted in a direction normal to 
laser beam. The spectra can be recorded either directly, 
or through the demodulation (with a Phase-Sensitive De- 
tection -PSD-) of a FM applied to the laser. Auxiliary 
laser beams allow the monitoring of the laser frequency: 
(i) one beam is sent to a scanning Fabry-Perot interfer- 
ometer to monitor (by means of PDA) the single-mode 
operation of the ECDL; (ii) the second one is used for an 
auxiliary saturated absorption (SA) set-up with a macro- 
scopic (3-cm long) Cs cell ensuring an accurate reference 
when scanning the ECDL frequency. 
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FIG. 1: Experimental set up. 



II. EXPERIMENTAL SET UP 



A scheme of the experimental setup is presented in 
FigfT] An extended cavity diode laser (ECDL) is used, 
performing frequency-tunable single-mode operation at 
A = 852nm, with a FWHM of about 3MHz. The main 
part of the laser beam, linearly polarized, is directed at 
normal incidence onto the ETC. The geometry of the ex- 
periment is chosen in a way that the laboratory magnetic 
field (about 0.5G) is approximately parallel to the laser 
light polarization. The construction of the ETC, filled 
with Cs vapor from a side arm, is similar to the one de- 
scribed in [4( . Its design was slightly modified to produce 
a wedge in the vapor gap. This makes the cell thick- 
ness locally variable in a convenient manner. The situa- 
tions L = A/2 or L = A are chosen by simply adjusting 
the relative position of the laser beam and of the ETC. 
The accuracy of the cell thickness measurement is better 
than 20nm. The Cs vapor density (~ 4.10 13 at. cm -3 ) 



III. EXPERIMENTAL RESULTS 
A. Sub-Doppler resonances in absorption 

The absorption spectra comprise two sets of h.f.s. com- 
ponents (a component being defined as optical transition 
between hyperfine sub-states) F g = 3 — > F e = {2,3,4} 
and F g = 4 — > F e = {3, 4, 5}. They are represented in all 
the following figures through the AP/Po ratio (denoted 
as absorption), where AP is the absorbed power, and 
Pq is the input power. The relative uncertainty on the 
transmission (Pq — AP) / Pq is on the order of a few 10~ 3 . 

Figure[5]illustrates the behavior of the absorption spec- 
trum on a cell of a thickness L = A/ 2, for three different 
irradiating intensities and for the F g = 3 — > F e = {2,3,4} 
set of transitions (a similar behavior is observed for the 
F g = 4 — > F e — {3,4,5} transitions). For all intensities, 
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the enhancement of the absorption at the center of the 
hyperfine transitions is responsible for a strong narrow- 
ing of the spectrum, notably allowing the resolution of 
the individual h.f.s. components. These results extend 
those of @, 0, [Hj]> evidencing the coherent Dicke nar- 
rowing; they are however obtained in an intensity range 
higher than the one (<< lmW/cm 2 ) ensuring a genuine 
linear behavior. At high intensities, the sub-Doppler res- 
onances appear significantly broadened, and saturation 
effects tend to washout the Dicke coherent narrowing [13] , 
which is well-pronounced at low power: one can notice in 
Figf2]that the absorption peaks are strongly reduced for 
high intensities, but that the wings are nearly unaffected. 
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FIG. 2: Absorption spectra across the F g 



{2, 3, 4} transitions for various intensities (as indicated) at 
L = A/2. 

Figure [3] shows the typical evolution of the absorption 
spectrum for L = A with the irradiating intensity (only 
the F g — 4 — > F e — {3,4,5} set is shown, but a similar 
behavior is observed for the F g — 3 — > F e — {2, 3, 4} tran- 
sitions). Significant differences are observed between the 
absorption spectrum at L = A/2 and that at L = A. Let 
us first recall that for L = A, no coherent Dicke narrow- 
ing is expected, and that in the linear regime, the absorp- 
tion profile, although complex, is Doppler-broad, owing 
to a (non velocity-selective) transient regime of interac- 
tion. Superimposed to the expected Doppler-broadened 
absorption profile, one observes, as a result of the rela- 
tively high intensities used here, well pronounced sub- 
Doppler narrow dips of reduced absorption. This re- 
duction of absorption is a signature of optical pumping 
and/or saturation processes that tend to reduce the num- 
ber of atoms available for the interaction with irradiating 
light. These processes can be completed only for atoms 
interacting a sufficient time with the laser light [l], 
they are highly enhanced for slow atoms (i.e. small ve- 
locity component along the normal to the ETC windows) , 
hence yielding sub-Doppler structures. The amplitude of 
these narrow structures increases when the irradiating in- 
tensity increases, the structures get apparently broader, 
and their contrast relative to the broad Doppler absorp- 
tion increases markedly, as the Doppler-broadened ab- 



sorption decreases under saturation. These results ap- 
pear very similar to those presented for the absorption 
spectra in ETC of Rb vapor [l2| , but include a regime of 
a higher irradiating intensity. 
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FIG. 3: Absorption spectra across the F g = 4 — > F L , = 
{3, 4, 5} transitions for various intensities (as indicated) at 
L — X. 
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FIG. 4: Experimental absorption at the line center of the 
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irradiating intensity for L = A/2 and L — A. 

Naturally, it is not a surprise that the narrow (velocity- 
selective) saturation dips are observed more easily when 
the non saturated lineshape is broad [L = A), than 
when it undergoes a notable coherent Dicke narrowing 
(L = A/2). To compare more quantitatively the differ- 
ing saturation behaviors for L = A/2 and L = A, the 
absorption at the center of individual h.f.s. component 
is plotted in Fig(3] as a function of the intensity. The 
comparison is here restricted to the F g = 4 — * F e = 4 
transition, but similar results are obtained for the other 
h.f.s. components, ft can be seen that the absorption 
rate decreases faster at L = A than at L = A/2. This 
faster reduction when the length increases could be seen 
as reminiscent of the behavior of velocity-selective pump- 



ing already observed in micrometric thin cells when 
the efficiency of the saturation process is governed by the 
product of the intensity by the cell length. An additional 
discussion is provided in Section IV. 



B. Narrow resonance in fluorescence 
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FIG. 5: Illustration of the ETC fluorescence spectra for lower 
(a, 7mW/cm 2 ) and higher (b, 130m W/cra 2 ) light intensities 
at the two studied cell thicknesses. Cs source temperature: 
127°C7 (a) and 105°C7 (b). 

As previously reported [!, 0], the fluorescence spec- 
tra exhibit sub-Doppler features that are narrower than 
those in the transmitted light, with an amplitude and 
width following a monotonic growth with the cell thick- 
ness. Very well resolved fluorescence spectra are recorded 
directly (Fig. O without FM and PSD of the signal. 
Even at low irradiation intensities (Fig. [5^,), the signals 
are narrower for L = A/2 than for L = A . At high irradi- 
ating intensities, we observe for L = A tiny dips that are 
superimposed to the top of the fluorescence profile of the 
individual h.f.s. components. Although these dips are 
observable through a direct detection (Fig. [5Jd), they are 
more conveniently characterized through the FM tech- 
nique (FigJSJ). Note as an additional difference between 
fluorescence and absorption spectra, that these satura- 
tion dips occur for much higher irradiating intensities (~ 
an order of magnitude in our experiments) in fluorescence 
than in absorption. 

To understand how these saturation features are spe- 
cific to fluorescence, two major differences with absorp- 
tion are worth being underlined : (i) while the absorption 
rate decreases to zero under saturation effects (because 
the population difference is reduced), increasing the ir- 
radiating intensity tends to increase the population of 
the excited state, and hence the fluorescence (at least 
as long the atomic system is a closed system); (ii) even 
in the linear regime (i.e. low irradiating intensity), the 
fluorescence is a second-order process [J] , that is velocity- 
selective in the transient regime specific to ETC. Hence, 
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FIG. 6: PSD fluorescence spectra (for L = A) as obtained 
after demodulation of an applied FM on (a) the set F g — 4 — > 
F e = {3, 4, 5}, recorded with an intensity : (1) -140mW / cm 2 ; 
(2) -1076mW/ cm 2 , and (b) on the set F g = 3 -> F e = {2, 3, 4} 
with an intensity : (1) -30QmW/cm 2 ; (2) -800mW/cm 2 . In 
the PSD spectra, the narrow dispersive structure exhibits a 
(quasi-) antisymmetry opposed to the one of the broad struc- 
ture. This is a signature of a narrow dip in the corresponding 
spectrum for direct detection. 



the observation of narrow dips in an already narrow flu- 
orescence spectrum can seem intriguing: in particular, 
for a closed atomic system, i.e. when no population loss 
occurs, the fluorescence rate is expected to be maximal 
at line center once the steady-state regime is reached. 
Moreover, one has to understand how a narrow velocity- 
selective dip of population loss can be superimposed to 
an already narrow velocity-selective peak of fluorescence. 
The last point can be tempered by the fact that the se- 
lectivity of atomic velocity decreases with increasing the 
cell length and is responsible for a significant broadening 
for lengths exceeding ~ A/2. This makes the width of a 
fluorescence spectrum not as narrow as the pure natural 
width and this leaves open the possibility of a more selec- 
tive process (population loss, assuming an open system) 
of the opposite sign; in addition, the sub-Doppler fluores- 
cence spectrum already undergoes a notable broadening 
because of the high intensities required to observe satu- 
ration dips. 

Before further interpreting our experimental findings 
(see Section IV) , it is interesting to point out that for the 
F g = 4 — > F e = {3, 4, 5} set of transitions, saturation dips 
appear only for the open transitions F g = 4 — > F e = 3 
and F g = 4 — > F e = 4, but are not observed for the 
closed transition F g = 4 — > F e — 5 (Fig. [6^,), in 
spite of the large explored range of irradiation intensi- 
ties (50 — lOOOmW / cm 2 ) . Conversely, Fig. [Sb strikingly 
shows that for the F g = 3 — > F e = {2,3,4} set of h.f.s. 
components, all h.f.s. components, including the closed 
transition F g = 3 — > F e — 2, exhibit a comparable satu- 
ration dip in the fluorescence spectrum under the consid- 
ered laser intensities. Actually, it is known that due to 
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the Zeeman degeneracy, a closed transition pumped with 
polarized light cannot be simply viewed as a transition 
on a (degenerate) 2-level system. It is in particular not 
protected against the Zeeman optical pumping that mod- 
ifies the tensorial orientation of the hyperfine sub-level 
through the excitation with linearly polarized light [3]. 
As a result of a strong irradiation on the F g = 3 — ► F e = 2 
transition, Cs atoms accumulate into the mp — ±3 Zee- 
man sub- levels, which do not interact with the laser light: 
although the transition is a closed one (in terms of energy 
level), a strong irradiation induces a decreased fluores- 
cence, because the system is actually an open one when 
considering the Zeeman degeneracy. At the opposite, on 
the F g = 4 — > F e = 5 transition, Cs atoms accumulates 
on Zeeman sub-levels with the largest absorption proba- 
bility mm, 

and a strong irradiation does not reduce 
the fluorescence. 



IV. DISCUSSION AND INTERPRETATION 
WITH A THEORETICAL MODELLING 

A. Limits in the interpretation and expectations 
from previous two-level models 

Although we describe a single-laser experiment, a fully 
quantitative prediction, would be very complex to obtain. 
This is because when dealing with saturation problems 
for degenerate two-level system, the Rabi frequency of 
an elementary transition between Zeeman components 
depends on geometrical Clebsh-Gordon type coefficients, 
so that a single parameter of saturation can hardly be 
defined. Moreover, here, the saturation process is gov- 
erned by a transient regime, and the duration of the in- 
teraction is velocity-dependent. This mixture of velocity- 
dependent transient regime, and of the velocity integra- 
tion, justifies that several regimes have already been an- 
alyzed in the elementary frame of a non degenerate two- 
level model. In all cases, the theoretical treatment of 
spectroscopy in a thin cell of dilute vapor assumes wall- 
to-wall atomic trajectory [E 0, S S G3 G3, El and the 
optical response results from a spatial integration of the 
local atomic response, that is determined through a tran- 
sient evolution. 

For relatively long cells, an elementary treatment of 
saturation, was developed P, Q (for an early indepen- 
dent approach, see [lq). relying basically on an open 
two-level model. Saturation effects were considered to 
be much slower than the coherent absorptive response, 
assumed to be instantaneous. They induce a velocity- 
selective dip in the absorption spectrum. An elemen- 
tary scaling law was found, with the key parameter (the 
pumping time) determined by the product " cell length by 
pumping intensity" . Even in this simplifying modelling, 
the identification of the atomic velocities contributing to 
the signal @ has revealed to be quite complex, because of 
an interplay between the velocity width associated to the 
natural optical width, and the maximal velocity allowing 



a quasi-steady state pumping. 

For lower intensities and/or smaller cell length, the 
optical pumping remains negligible, and an elementary 
(closed two-level) model has to be considered @. The 
relevant transient regime is the build-up of the absorb- 
ing properties of the vapor, i.e. of the optical coherence. 
The interference between the various velocity-dependent 
(complex) coherent response of all atomic velocities is at 
the origin of the Dicke narrowing at L = A/2 (and of 
its periodical revivals at L = (2n + l)A/2). The Bloch 
vector model is an adequate tool to explain the pe- 
riodical Dicke narrowing for absorption and it can be 
applied beyond the limits of the linear regime [H, l2C| - 
or to accommodate the population losses of an open sys- 
tem. For a strongly driven irradiation, and a moderate 
relaxation (closed system) , the global process remains 
purely coherent, but requires the velocity integration (or 
interference) of quickly rotating Bloch vectors. Through 
the interference of these multiple oscillations, it can be 
inferred [H, H3| that varying the detuning, (i.e. chang- 
ing the orientation of the pseudo-magnetic field in the 
Bloch- vector model), will lead to an oscillating behavior, 
instead of always yielding a maximum at line center. And 
indeed, a formal analytical treatment (Toj | for a two-level 
model (closed or open system), predicts such a multi- 
peaked absorption lineshape under a strong saturation, 
starting with a simple dip at line center for moderate sat- 
uration. However, in spite of its formal analyticity, this 
treatment requires a numerical determination of the rele- 
vant eigenvalues determining the solutions, and becomes 
cumbersome for a lineshape calculation. Also, no specific 
analysis has been provided for the situations the most rel- 
evant for the coherent Dicke narrowing, only a situation 
closed to L — 3A/2 is explicitly studied. In particular, for 
an open system, the problem of the competing physics 
of the coherent response in the saturation regime, and 
of the incoherent velocity-selective population depletion, 
has not been addressed in [lOj. The systematic mod- 
elling of the fluorescence response in the context of ETC 
has never been reported. In [131 ], fluorescence spectra in 
ETCs are calculated, but they mostly aim at the specific 
description of the multi-level Rb transition. 

As a further step of our study, we propose a numer- 
ical evaluation of the response of an elementary non- 
degenerate two-level atomic system (conservative or open 
2-level system), and we compare it with the experimental 
findings. 



B. A two-level modelling in view of a numerical 
estimate 

We consider here a two-level model with a control pa- 
rameter allowing to compare the situation of an open 
system (with losses to a generic "third level" 1'), or of a 
closed system (FigJT]). Levels 1 and 2 are coupled by a 
laser light at a frequency u, detuned by A from the tran- 
sition frequency u>2i (A = u>— u>2i)- The Rabi frequency is 
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defined as = 1\i\zE§lYi {Eq being the light field input 
amplitude) with /i 12 the dipole moment of the transition. 
The width of the excited state is denoted by 72 , and one 
assumes, for the considered dilute vapor, that collisions 
- notably dephasing collisions- can be neglected, so the 
optical width of the transition 721 is 721 = 72/2. Prac- 
tically, this simple hypothesis introduces a set of dual 
relaxation constants. This unfortunately leads to hardly 
tractable analytical solutions, while the assumption of a 
single relaxation constant for population and optical co- 
herence, would have greatly simplified the calculations. 
However, such a simplifying assumption can hardly be 
justified in the context of an ETC, with the wall-to-wall 
atomic trajectories 21]. To take into account the pos- 
sibility of population losses to a third level (e.g. to the 
other hyperfine sub-level of the ground state of alkali- 
metal atoms, or to the Zeeman sub- levels), one intro- 
duces a coefficient, defined as a, for characterizing the 
probability to decay from level 2 to level 1 (a = 1 for a 
closed system, < a < 1 for an open system). 
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processes that are altogether essential for the sub-micron 
sized cells. 

As usual in a thin cell, the initial conditions for 
the system of Eqs([T][3]) differ for arriving (v < 0) 
and departing (v > 0) atoms. One has indeed 
[cr n (L) = l;a 22 (L) = 0;o 2 i(£) = 0] for v < 0; and 
[on (0) = 1; 0-22(0) = 0;o2i(0) = 0] for v > 0. To relate 
the solution of the system (1-3) to the signals of ab- 
sorption or fluorescence observed in the experiments, we 
further follow the approach presented in Q. The local 
atomic response (at z) is deduced from the integration of 
its transient behavior (owing to z = vt, or z = L + vt, for 
respectively v > and v < 0). The optical signal results 
from the spatial integration of the atomic response, af- 
ter the required integration over the velocity distribution 
(assumed to be a Maxwellian, with a thermal velocity u). 
Hence, the absorption is proportional to a quantity 
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G{v)exp[-{^f]dz 

[ j KVb 



with 



G(v) 



Im[u2i {z, v)]dz. 



(4) 



(5) 



In the experiment, the measured signal is the coher- 
ent beating between the input field and the reemitted 
field It ~ 2EoEt, so that it is the experimental ratio of 
absorption AP/Pq which has to be compared with the 
theoretical quantity A/ flu (or It/Po). In a similar way, 
the fluorescence is related to the quantity U, with : 



FIG. 7: Schematics of the atomic system used in the theo- 
retical modelling. 

The system of Bloch equations is hence the following: 



1=1 Q(v)exp{-A 2 }dv 
Iq ku 



(6) 



where Q(v) is defined as 



v— I--D21C21 - i-^- (011 - 022) = 0, (1) 

az 2 



Q(v) 



[0-22(2, v)]dz. 



(7) 



d<7> 



22 



dz 



dan 

dz 



72CT22 - Q.Rlma 2 i = 0, 



072022 + Vl R Ima 2 i = 0, 



(2) 



(3) 



where D 2 i = 721 + ikv — iA , v is the atomic velocity 
(along the laser beam, and hence along the normal to 
the ETC), and the aij are the reduced density matrix el- 
ements in the rotating frame. The above system had been 
solved analytically [l|, in an approach focusing only on 
the non-linear incoherent processes. This could be jus- 
tified in view of solving the restricted problem for cell 
lengths allowing transient coherent processes to be neg- 
ligible. Here we consider both coherent and incoherent 



On this theoretical basis, it is possible to spatially- 
and velocity- integrate the solutions of the density ma- 
trix equations to be found numerically under the condi- 
tions of saturation. Relatively to the previous calcula- 
tions PISH) one should recall that: (i) the modelling of 
the coherent Dicke narrowing was achieved on the basis 
of similar density matrix equations 0, [13] , in the limit of 
a first-order interaction with the resonant light (in such 
a case, the parameter a, affecting population redistri- 
bution, plays no role in this first order prediction, and 
the Bloch vector model applies), and that (ii) for the 
velocity-selective optical pumping in thin cells (L >> X 
) [l|, the optical coherence yielding the absorption rate 
was estimated under a rate equation approach, allowing 
the instantaneous measurement of the remaining active 
population difference (for a < 1). 
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In view of discussing some of the theoretical predic- 
tions with parameters applicable to an elementary and 
realistic case, the above model has been used with the 
following parameters: 721 = 5MHz, ku = 250MHz, 
and a — 1 (closed system) 01 a — 0.5 for a realistic open 
alkali system. Technically, our numerical results combine 
a velocity integration, the spatial integration of a locally- 
varying response (see Eq[5]and Eq[7]), and a Runge-Kutta 
integration equivalent to the integration of the transient 
response governing the spatial response for a given veloc- 
ity. 



Comparison between the model and the 
experiments : Absorption behavior 
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FIG. 8: Intensity dependence of the absorption - at line- 
center- for a closed transition. The theoretical dependence 
for L — A/2 is in dashed line, the solid line is for L = A. 
The experimental data {L = A/2 : open squares, L = A : full 
squares) are for the closed F g = 4 — > F e = 5 transition. For 
the scale applied to the theoretical curves, see text. 



nents because they are observable only over the slope of 
a stronger component, adding an extra-difficulty to char- 
acterize the appearance of an inverted dip structure. A 
rigorous measurement for a given hyperfine component 
would imply to subtract the contribution from the neigh- 
boring components. Such an evaluation cannot be very 
precise, and it becomes natural to concentrate the study 
on the stronger transitions. One can also mention that 
the uncorrected transverse structure of the irradiating in- 
tensity (presumably Gaussian) tends to wash out the tiny 
oscillations that could be induced by saturation [Tfj| , and 
that the residual terrestrial magnetic field, although not 
sufficient to generate a resolved Zeeman structure, may 
modify the coupling rates between sub- levels. 
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When attempting to compare quantitatively the ex- 
perimental results with a modelling, and especially if it 
is not intended to go to a complete lineshape analysis, it 
is necessary to recall various intrinsic limitations affect- 
ing the possibility of a quantitative comparison between 
the experiment and the above two- level model. First, 
this comparison remains in principle of a limited scope 
because multiple Zeeman transitions are involved. This 
means that in principle, saturation effects cannot be ac- 
counted by converting the experimental intensity into a 
single Rabi frequency Qr. The tensorial structure of the 
atomic system makes non identical the various transfer 
rates to the individual sub-levels. On the experimental 
side, the hyperfine components are not perfectly resolved, 
but partially overlap, moreover in a non constant man- 
ner that depends upon the cell length, and the saturation. 
This makes uneasy to attribute all of the measured ab- 
sorption at a given frequency - on the center of a h.f.s. 
component, or elsewhere - to a single h.f.s. component. 
This limitation is even stronger for the smallest compo- 



FIG. 9: Comparison between the intensity dependence of 
the absorption on a open transition, at line-center, and at 
a slightly detuned frequency. The ETC thickness is (a) 
L — A/2; (b) L = A . The theoretical dependence is shown in 
solid line for a frequency at line center A = 0, and in dashed 
line for a detuned frequency A = O.OSku. The experimental 
data (A = +20MHz : open squares, A = : full squares) 
are for the open F g — 4 — * F e — 4 transition. For the scale 
applied to the theoretical curves, see text. 

For all the above reasons, the interplay between nu- 
merous processes, with differing time constants, makes 
hopeless the characterization of the complex broadening 
of lineshapes by a "width" of resonance. This is why, 
in an attempt to simply evaluate the onset of the ap- 
pearance of a narrow inverted structure (i.e. reduced ab- 
sorption), we compare predictions for an irradiation fre- 
quency at line-center, and for a slightly shifted frequency 
(we take A = 0.08fcu, or 20MHz for numerical values 
as mentioned in Section IV-B). Such a criterion, possi- 
bly misleading if the spectrum would include numerous 
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oscillations, seems reasonable with respect to the appar- 
ent width of the various saturation dips that we observe. 
Figsl8lfT0l allow a comparison between the experimental 
and the theoretical results for the F g = 4 — > F e = 4 
and F g = 4 — > F e = 5 transitions. These two transi- 
tions are good examples of open and closed transitions, 
and the corresponding "saturating intensity" (although 
the concept is, as mentioned, of a limited scope for a de- 
generate system) should be quite comparable. To make 
the theoretical predictions directly comparable to the ex- 
periments, we use a conversion factor (£Ir/j2i) 2 = 1 for 
15mW/cm 2 which was chosen to provide the most satis- 
factory visual fit between the experimental and theoreti- 
cal curves. Also, in these figures, the absorption rate for 
the theoretical curves was adjusted (by a factor of 1.8) 
to provide the optimal comparison with the experiments; 
note however that in principle, the absorption rate is pre- 
dictable in an absolute manner provided that the atomic 
density and the dipole moment are known. 

This comparison between the simplified modelling and 
the experimental observations shows a satisfactory agree- 
ment. In particular, in Fig. [8) where the predicted ab- 
sorption at line center is plotted in the two typical cases 
L = A/2, and L = A , one notes as predicted that if the 
absorption at L = A/2 is only slightly smaller than for 
L = A for low intensities (as expected due to the Dicke 
narrowing, the exact ratio being governed by the 721/fcu 
factor), the absorption becomes even larger for L = A/2 
than for L = A at higher intensities. Interestingly, such 
a result is valid for a closed system (Fig. [8]) as well as 
for an open system - see FigJHl Although the saturation 
processes for L = A that reduce absorption of slow atoms 
are in principle twice more efficient than the compara- 
ble processes for L = A/2 , the dominant effect seems 
here to be the survival of the coherent Dicke narrowing 
(for L — A/2 ), with its large contribution of fast atoms 
that are nearly insensitive to the saturation. This larger 
contribution at line center is the distinctive evidence of 
the coherent Dicke narrowing, induced by the coherent 
transient contribution of atoms that are not " slow" . It is 
hence natural that the Dicke coherent narrowing remains 
quite robust, as unaffected by relatively strong irradiat- 
ing intensities. However, it cannot be concluded that a 
narrow saturation dip in the absorption would not be 
observed in the conditions allowing for a revival of the 
Dicke narrowing, such as L = 3 A/2 (a length unfortu- 
nately not attainable because of the construction of our 
cell): indeed, in most experimental conditions, the re- 
vival of the Dicke narrowing [|| @] (although shown to be 
robust with saturation at L = 3 A/2, see pjl) only brings 
a sub-Doppler structure of a small amplitude. 

Figures 151 and ITD1 allow the comparison of the L = A/2 
and L = A situations, with respect to the appearance of 
a narrow dip of reduced absorption, for closed and open 
systems. For the closed system considered in FigflOl one 
predicts (FigllOb) for L = A that the initially broad peak 
- i.e. no Dicke narrowing- exhibits an inverted substruc- 
ture even for low intensities (a few mW/cm 2 ) . Con- 
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FIG. 10: Same as Fig. but for a closed transition theory 
(experimental data from the F g = 4 — » F e — 5 transition). 

versely, for L = A/2 (FigllOb,). the narrow Dicke struc- 
ture undergoes only a visible broadening, but without the 
clear appearance of a dip in the center of the transition. 
A closer look on the inset of Fig. [TUb shows however 
that A = is no longer the peak of absorption for high 
intensity, but the amplitude of the corresponding dip is 
predicted to be extremely small. This demonstrates that 
the absence of observation of a narrow dip for L = A/2 
is not some fundamental effect, but rather the quantita- 
tive result of the competition between distinct processes 
affecting optical coherences (Dicke narrowing) or atomic 
population (saturation). An analogous behavior is pre- 
dicted for an open system (FigJSJ a = 0.5), with sat- 
uration at line center and a tiny dip for L = A/2 at 
high intensities (slightly more pronounced than for the 
closed system) , and the occurrence of a pronounced nar- 
row dip for L = A. For L = A, the dip amplitude is 
predicted to be significantly larger for the open transi- 
tion than for the closed one which is not observed in the 
experiment. For these discrepancies, it should however 
be kept in mind that our analysis here tackles narrow 
details of the lineshapes, and that a shift of 2QMHz is 
experimentally small. 

D. Predictions for the fluorescence behavior 

Because saturation effects in fluorescence are observed 
for higher intensities than in absorption, and mostly in 
a PSD technique following an applied FM, we have not 



9 



attempted to perform a quantitative comparison between 
the experiments, and the predictions of the modelling. 
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FIG. 11: Calculated fluorescence (1,2) and absorption (3,4) 
profiles for open (1,3) and closed (2,4) transitions under condi- 
tions where reduced absorption dip is predicted: = 8. 

The numerical calculations confirm the experimental 
observation that lineshapes are expected to be narrower 
in fluorescence than in absorption in comparable condi- 
tions [J, lZ|, |l3j. Also, the width of the fluorescence profile 
is expected to increase continuously with the cell thick- 
ness, without an interferometric Dicke-type narrowing, 
and to reach a Doppler-broadencd lincshapc for longer 
cells with the velocity selection getting less stringent. 
The simulation in FigQj] (i.e. calculation with the rel- 
evant experimental parameters) does not predict the for- 
mation of a narrow dip in fluorescence for L — A/2 and 
this agrees with our experiment (Fig. [5]), while under 
the same condition, a dip in the absorption is predicted. 
If such an absence of a dip can be expected for a closed 
system, a strong irradiation should be able to induce a 
severe depletion of the fluorescent atoms for an open sys- 
tem, and hence a dip in the lineshape. However, in the 
sense of the dip formation, our simulation does not show 
essential differences between the closed and open tran- 
sitions (Fig JIT]). Most probably this is because, for our 
choice of 712 / ku parameter this strong irradiation would 
imply for L = A/2 such a large broadening of the tran- 
sition that the velocity-selection itself is governed by a 
width not markedly narrower than the one of the total 
signal. 

When increasing the cell length to L = A, Fig[T2"a,b 
shows, with a presentation similar to the one used in 
Figs I9TTUI that a pronounced dip of reduced fluores- 
cence at line center is predicted at L = A for the open 
transition. For a closed transition, no dip is predicted, 
and rather, the lineshape broadens with saturation. In 
addition, theoretical fluorescence profiles are presented 
in Fig|T2"b, showing completely different behavior of the 
fluorescence around the transition center for the open 
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FIG. 12: Theoretical features of fluorescence for L — X. Inten- 
sity dependence of the fluorescence (normalized by the Rabi 
frequency, i.e. plot of U ) at line center A = (solid line), 
and at a slightly detuned frequency A = 0.08fcit (dashed line) 
for : (a) an open transition; and (b) a closed transition, (c) 
Theoretical lineshapes for open (1) and closed (2) transitions 
for fifl/721 = 2. 



and closed transitions. This striking theoretical differ- 
ence between closed and open system justifies our exper- 
imental observations (Sect. III-B, and FigO, where a 
fluorescence dip is observed for all components but the 



F„ 



5 transition. 



V. CONCLUSIONS 

In spite of the apparent simplicity of single beam ex- 
periments on ETCs, we demonstrate a large variety of 
regimes in the study of saturation. This is because 
thin cell spectroscopy naturally yields a signal averaged 
on various regimes of transient interaction between an 
atomic velocity group, and a resonant irradiation. This 
also explains that a large variety of modelings has been 
proposed to deal with these effects. As long as the real 
system of an alkali atom (such as Cs in our case) is far 
from being a two-level system, owing to its degenerate 
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multi-level nature and including the tensorial structure 
responsible for the various Zeeman sub-states, it is hope- 
less to describe in full detail the saturation effects : this 
can be easily understood by recalling that for alkali va- 
por, there exists no general description of saturated ab- 
sorption spectra under a strong pump irradiation: this 
latter problem is however notably simpler as being lim- 
ited to a steady state interaction, but similarly sensible 
to the many coupling strengths involved in the highly 
complex sub-Doppler atomic structure of alkali atoms. 
In this context, it becomes clear that an exact quantita- 
tive description would be an enormous task of a proba- 
bly limited benefit. Nevertheless, it is remarkable that 
a comparison between a pure two-level model, with well- 
chosen numerical parameters, and our experiments, leads 
to a relatively satisfactory agreement. 

On more general grounds, some major features can be 
deduced from our studies. In absorption, the satura- 
tion reduces preferentially the contribution of slow atoms. 
This leads to the observation of saturation dips on line 
centers, that are observed more easily when the (non sat- 
urated) absorption is broad (i.e. for L = A ), than when 
the coherent Dicke narrowing makes the lineshape intrin- 
sically narrow. For L = A/ 2, the Dicke narrowing is so 
robust that we only observe a broadening, without the oc- 
currence of the predicted tiny saturation dip. Because of 
the coherent nature of absorption processes in ETCs, the 
observed dips at line centers can result from the combined 
velocity-selective depletion of population difference, and 
from a complex oscillating behavior. These oscillations 
are probably more efficient for ETC thickness leading to 
a coherent narrow Dicke structure, than for lengths mul- 
tiple of L = A, characterized by destructive interferences 
across the Doppler-broad structure, and for closed sys- 
tems rather than for open systems (as characterized by 
an incoherent velocity-selective population transfer). 

In fluorescence, it is only for open atomic systems, al- 
lowing a reduced contribution of slow atoms, that dips 
at line center can be observed as originating from a 
velocity-selective process. Also, increasing the cell length 
makes easier the observation of a narrow structure in- 
side the sub-Doppler fluorescence spectrum as due to the 



increased width of the non-saturated fluorescence spec- 
trum. For effective closed systems (i.e. including the re- 
distribution among the Zeeman sub-levels), it is not clear 
if the fluorescence spectrum can involve oscillations, rem- 
iniscent of the kind of fringes that are predicted to ap- 
pear in the absorption spectrum. This possibility could 
strongly depend upon the relative relaxation of the opti- 
cal coherence rate, and population losses. In our exper- 
iments, the relatively strong coherence losses, owing to 
the relatively high Cs temperature, and the uncontrolled 
spatial distribution of the irradiating beam, could be suf- 
ficient reasons to make unobservable saturation features 
more complex than an elementary dip. 

The reported results enhance our knowledge in the 
rich field of the Doppler-free ETC spectroscopy which 
is of significant importance for the development of high- 
resolution spectroscopy of atoms and molecules confined 
in nano-volumes. ETC spectroscopy has recently been 
shown to allow the spatial analysis of the long-range van 
der Waals atom-surface attraction, that modifies spec- 
tra for short ETC thicknesses (22J. The recent obser- 
vation [23j of Electromagnetically Induced Transparency 
effect in ETC is promising for the dynamics study of 
this widel y us ed phenomenon. ETC application has been 
proposed [2J] for magnetic field measurements with sub- 
micrometer spatial resolution which can be useful for de- 
tailed magnetic mapping performance. 
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